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Abstract 

A variety of infiltration techniques can be used to fabricate solid materials, particularly 
composites. In general these processes can be described with at least one time dependent 
partial differential equation describing the evolution of the solid phase, coupled to one or 
more partial differential equations describing mass transport through a porous structure. 
This paper presents a detailed mathematical analysis of a relatively simple set of equations 
which is used to describe chemical vapor infiltration. The results demonstrate that the 
process is controlled by only two parameters, a and 0 . The optimization problem associated 
with minimizing the infiltration time is also considered. Allowing a and 0 to vary with time 
leads to significant reductions in the infiltration time, compared with the conventional case 
where a and 0 are treated as constants. 
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1 Introduction 


A variety of materials are produced by infiltration processes. In these techniques a fluid 
phase (i.e., a gas or a liquid) is transported into a porous structure, where it then reacts to 
form a solid product. These methods are particularly important for producing composite 
materials, where the initial porous preform is composed of the reinforcement phase (i.e., 
fibers, whiskers, or particles) and infiltration produces the matrix. ([1]),([2]). A detailed 
assessment of the relevant reaction and mass transport rates during infiltration requires 
mathematical modeling, using a minimum of two coupled partial differential equations which 
describe changes in the reactant concentration and the solid structure as a function of both 
position and time. This type of modeling can also be extended to analyze the optimization 
and control of infiltration processes. 

The research presented here specifically considers optimization for a set of two equations 
which describe chemical vapor infiltration (CVI). In this process a vapor-phase precursor is 
transported into the porous preform, and a combination of gas and surface reactions leads 
to the deposition of the solid matrix phase. In recent years a number of researchers have 
developed mathematical descriptions of CVI. While these models do not provide complete, 
detailed representations of CVI, they provide an excellent starting point for mathematical 
and computational research on the equations which describe infiltration processes in general. 
During infiltration the formation of the solid product phase eventually closes off porosity at 
the external surface of the body, blocking the flow of reactants and effectively ending the 
process. This is a key feature of most infiltration processes. It is usually desirable to 
maximize the amount of solid formation that occurs before this endpoint is reached. This 
optimization goal is often linked to another objective. For example, CVI often requires 
extremely long times, so it is also important to minimize the total processing times. 

The paper is organized as follows: In Section 2 we derive a simple set of two equations 
that models the infiltration processes, and we show how to get these equations as a subcase 
of the general system. These are nonlinear partial differential equations. We also derive 
initial and boundary conditions for the model equations. In Section 3 we present a detailed 
mathematical analysis concerning the behavior of the solutions in space and time. This 
analysis provides insights into the behavior of the process. In Section 4 we define the concept 
of a successful process and get conditions on the paremeters of the problem for a process to 
be successful. In particular we formulate an optimization problem for the minimum time it 
takes for the process to settle. In Section 5 we present numerical experiments to validate the 
theory. We also discuss how to design the experimental parameters to get faster successful 
processes. 


2 Formulation 

A mathematical description of infiltration requires one or more partial differential equa- 
tions which describe the evolution of the matrix (i.e., the solid phase), and one additional 
partial differential equation for each chemical species in the fluid phase. For a simple pore 
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structure, the continuity equation for species i is 

?<|3i + V • JV, = £>, r -R, (1) 

at r 

where t is time, £ is the void fraction of the media, C t is the concentration of species i, n r 
is the number of the gaseous species, u ir is the stoichiometric coefficients for the ith gaseous 
species in the rth reaction, and Rr represents the volumetric reaction rate of reaction r. 

The basic partial differential equation (s) which describe reaction and mass-transport 
in porous media (i.e., the fluid phase) are well-established [6, 8]. The Dusty-Gas model 
[9] describes multicomponent diffusion and convection. Neglecting thermal diffusion, the 
relationship between the molar fluxes, N t , is given by [10]: 


Nj RT ^ C jNj CjNj _ CjBe ^ p 

D Ki P V Mi} ' fiD Ki 


( 2 ) 


where B e is the permeability of the porous media, // is the viscosity of the mixture, and P is 
the total pressure. V MiJ and D Ki are the effective binary diffusivity for species i and j and 
the effective Knudsen diffusivity of species i, respectively. 

The change in the solid structure is equivalent to considering the change in the void 
fraction, e (i.e., the volume fraction of gas inside of the porous solid). The evolution of £ is 
given by: 


= —uS v (e) (3) 

at 

where u is the rate at which the solid product grows (volume/area/time) and S v (e) is the 
gas/solid surface area per unit volume of the porous solid. 

The simplest formulation for the fluid phase is obtained by considering one reacting 
species. For highly diluted reactant systems, the Dusty-Gas model can be simplified to give 
the following approximate expression for the flux: 


jv =-°i (4 » 

where C is the concentration of diluted species and Z is the distance into the preform. The 
effective diffusivity of the diluted species, D, can be expressed as 


D — - Dm m [1 + A'fc(c)] 1 (5) 

where m refers to the bulk species, D Mm is the binary diffusion coefficient for M in m, 
N k is the ratio of the Knudsen diffusion coefficient and and 0 is the tortuosity factor. 

For a diffusion-limited process in one spatial dimension, using Eq. (4), Eq. (1) becomes: 


d{sCi) _ d [ dCl uS v (s) 

dt dZ [ dZ J V M 

where V M is the molar volume of the solid product. The last term in Eq. (6) describes the 
rate at which the gas-phase precursor is consumed (or created) by chemical reactions inside 
of the pores, with the assumption that there are no homogeneous gas-phase reactions. 
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A specific CVI model requires expressions for u, S v , and D. Our objective in this work is 
to use simple formulations for each, as a basis for assessing the general behavior of infiltration 
problems. As an example, consider the formation of carbon matrix composites using a 
hydrocarbon in an H 2 carrier gas, where the following net reaction occurs: 

C m H„ (g) ^ m C(s) + | n H 2 (g) (7) 

The form of Eq. (6) is based on the assumption that the MTS concentration, C r , is dilute 
(i.e., the reactant concentration is much smaller than the carrier gas concentration). If the 
carbon growth rate is proportional to the precursor concentration, then: 


u = kC r (8) 

where k is the reaction rate constant. 

The preforms used for CVI typically have a complex porous structure. However, a cylin- 
drical pore is often used to formulate simple models. This leads to the following expression 
for S v : 


Sv{e) = 


2 yj£~o\f£ 


where r 0 is the initial pore radius and e 0 is the initial concentration of e. 

Substituting Eqs. (8) and (9) into Eqs. (3) and (6) gives the following forms: 


d V 

m = ~i Bc 


d_ 

dz 


[1 + AM)]' 


.1 dc 
dz 


where: 


a 2 r)c 


( 9 ) 


( 10 ) 

( 11 ) 


t] - sfe 

(12) 

C r 


c ~ cl 

(13) 

z 

(14) 

Z ~L 

q2 _ 2 kjr 0 L 2 

V M ToD M m 

(15) 

C lkyfe 0 C 0 2^ mC'oDmTTI 

r 0 ~ a L 2 

(16) 


where L is the half-thickness of the preform, C Q is the concentration of the reactant species 
in the bulk gas-phase (i.e., outside of the preform). The expression for a (15) is based on 
the assumption that u is determined by a first order reaction, where k, is the rate constant 
(i.e., u = kC r ). Note that a 2 is dimensionless and that j3 has units of inverse time. 
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The parameters a 2 and 0 depend on the three key process variables: T, P , and C 0 . T, P 
do not appear explicitly in Eqs. 15 and 16, however, k typically obeys an Arrhenius-type 
exponential temperature dependence, and D Mm varies with both temperature and pressure. 
If the process variables are all held at single fixed values throughout the process, then a 2 
and 0 differ only by a lumped constant. However if one or more of the process variables is 
changed as infiltration proceeds, then the proportionality between a 2 and 0 wdll also vary. 
As an example of this, we consider a case where a 2 decreases while 0 is held constant (see 
Section 5). One wav of accomplishing this in practice is to reduce the temperature during 
the course of the process, such that a 2 decreases (because k usually decreases faster than 
D M m as T decreases). To hold 0 constant, C 0 must be increased in a way which directly 
offsets the decrease in k. 

The time derivative in Eq. (6) has been removed in Eq. (11). This is permissable because 
solids are much denser than gases, such that the time-scale for changes in the C T profile is 
much shorter than the time scale associated with changes in the solid structure. For gas- 
solid reaction processes such as CVI, this is sometimes referred to as the pseudo steady-state 
approximation [7]. Transforming £ to T] simplifies equation (10). Basically tj is propotional 
to S v , so it is also possible to view 7? as a dimensionless surface area per volume. 

The boundary conditions that are most often used for CVI models are to fix the con- 
centration at the outer surface of the preform at C 0 , and to assume that diffusion occurs in 
from two opposite sides, such that there is no net flux in the middle of the preform (i.e., at 
Z = L): 


c(0,i) = l (17) 

c*(l,t) = 0 (18) 

The initial condition is given by: 

e{z.S)) = e 0 (19) 

During CVI, the infiltration kinetics are controlled by diffusion and the deposition reac- 
tion. To achieve relatively uniform infiltration, diffusion must be fast relative to the deposi- 
tion rate. This is typically accomplished by choosing processing conditions that result in a 
slow deposition rate, which usually leads to long infiltration times. Thus, a key processing 
objective is to obtain the desired amount of infiltration in the shortest possible time. The 
total amount of infiltration in the preform is given by integrating over 2 : 

e{t) = f e{z,t)dz (20) 

J O 

It is generally important to obtain the desired density (i.e., void fraction), £f, in the 
shortest possible time. Thus, the optimization problem of interest corresponds to determin- 
ing the shortest time where e(t) = £/, for values of e f that are significantly smaller than 

So- 

3 Properties 

The following system is obtained from Eqs. (10) and (11): 
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( 21 ) 

( 22 ) 


_ -Ml c(z t) 

at 2 1 ’ ; 

= a 2 (f)?7(2,«)c(^t) 

where f(r]) = \ [1 + iW^)] -1 , € C°° is a monotone increasing function of 77 , a(t) ^ 0 

and 3(t) > 0 are C°° functions of t, and 77 is yfi. 

(23) 

(24) 

(25) 


The system is subject to initial conditions 

77 ( 2 , 0 ) = 7/0 > 0 

as well as boundary conditions 


c( 0 . t ) = 1 


dc 

dz 


(1 ,t) = 0 


Physically, these correspond to a constant concentration at the outer surface (2 = 0 ) and 
a symmetry plane at 2 = 1 . We can immediately find the initial value for the concentration 
function c(z, t ) 


Lemma 1 

The initial distribution of the function c(z, t) is given by: 


where 


c(z, 0) 


cosh 7 o(l — 2 ) 
cosh 70 


7o = 



(26) 


(27) 


Proof 

Consider ( 22 ) at the initial time t = 0. Because of the initial condition (23) 77 ( 2 , 0) = rj 0 is 
a constant in 2 . Equation ( 22 ) is, therefore, a second order constant coefficient ODE whose 
solution is given by (26). 


□ 


Lemma 2 

The value of 77 ( 2 , t ) at the boundary point 2 = 0 is given by 

77(0, t) = 3 o ~ ^ P(r)dr (28) 


Proof 

We read ( 21 ) for 2 = 0 using the boundary condition (24) c( 0 .f) = 1 . Integration yields 
(28). 
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□ 

Examining the boundary value for 77(0, t) given in (28). it is clear that a positive solution 
does not exist if 

\ [ P(r)dT = Tjo (29) 

2 Jo 

for some t. We define the first time that (29) is satisfied as the terminal time tp. If (29) is 
never satisfied we define the terminal time tp as 00. 

We assume now that before the terminal time, 77(2 , t) exists and is positive, to be precise 

Assumption 1 

There exists a unique positive solution 77(2, £) for any 0 < t < tp, 0 < z < 1. This 
solution is C°° in both 2 and t. 

□ 


Comment: 

7, From Assumption 1, it follows that there exists a solution c(z, t) in the same time and 
space intervals. 

In the rest of this section we will derive properties of the solutions c(z,t) and 77(2, f). 
Study of these properties is important not only for the understanding of the behavior of the 
infiltration process, but also to find ways to speed up the process, in particular it will help 
us to define a successful process. 

In the following Lemma we will consider the behavior of the concentration function c(z , t ) 

Lemma 3 

The concentration function c(z, t ) is a positive function for 0 < 2 < 1 for any 0 <t <tp. 
Moreover, c(z, t) decreases as a function of 2. 

Proof 

We multiply (22) by c(z , t ) for every 0 <t<tp and integrate from some point 2 to the 
point 2 = 1 taking into account that the spatial derivative of c(z, t) vanishes at 2 = 1. For 
convenience we drop the dependence on t. We then get 

-c W /(„ W )^ = /'/««) 

+<x 2 J z > ( 30 ) 

thus c{z) is negative. It can never vanish, since, if c(z) vanishes at a point z 0 then 
the right hand side of (30) has to vanish and thus both c and ^ vanish in the interval 
[2 0 , lj. Bv the standard theory of linear differential equations the solution c(z, t) must vanish 
everywhere. So that c(z)^~ is strictly negative. Since c(0, t) = 1, it follows that c(z ) is 
positive and is negative. Thus the lemma is proven. 
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□ 

We turn now to the behavior of rj(z, t) as a function of the spatial variable z, for any 
fixed t. 

Lemma 4 

The function rj(z, t) increases as a function of z. In particular 

V(0,t) < rj(z, t) < 770 (31) 

where 77(0, f) is given in (28). 

Proof 

Consider equation (21) and integrate with respect to t to get 


1 f l 

r)(z, t ) = 770 - - J 0(t)c(z, r)dr (32) 

By Lemma 3, c(z,t ) is monotonically decreasing (in z), therefore rj(z,t) is monotonically 
incre«ising (in z). 

The lower bound in (31) is a consequence of the monotonicity of 77, the upper bound is 
a result of (32) and the fact that c(z, t) is positive. 


We can also establish an upper bound on the concentration function c(z,t ): 

Lemma 5 

The concentration function c(z, t ) is bounded by 


where 


c(z,t) < 


cosh 71 ( 1 — z) 
cosh(7!) 


□ 


(33) 


71 = \f^i ait) (34) 

In particular if a (t) = q( 0), is independent of time, then 

c(z,t) < c(z,0) (35) 


Proof 

Denote first bv F(z, t ) the solution of 

J^ F (z,t) -l\F{z,t) = 0 


(36) 
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with the boundary conditions 


F(0,t) = 1 


lr <M) = 0 ' 


The solution of (36) is clearly 


F(z,t) = 


cosh 7i(l — z) 
cosh(7x) 


(37) 


(38) 


We rewrite now equation (22) as 


d 2 c(z,t ) f z dc(z,t ) 

dz 2 + / az 


7 2 c(z,t) 


(39) 


where 


7 = 


N 7fa(M)) 


a{t) 


(40) 


We multiply now (36) by c(z,t), (39) by F{z,t) and integrate between some z and 1, 
taking into account that the spatial derivatives of both c(z,t) and F(z,t) vanish at z = 1, 
to get: 


~( F §-f) = -fyt F W> 2 -^ (41) 

Now f z > 0 since / is monotone in tj and r] is monotone in z. By Lemma 3, ^7^ is 
negative and so while F is positive. Also we have assumed that ^ is monotone so that 
7 2 — 7 2 is positive. We conclude that the RHS of (41) is positive. Thus 


Be dF 
F-tt ~ c-z- 
dz uz 


< 0 


(42) 


Leading to 


d c(z.t) 

F[ZT) 


Since at z = 0. c(z,f) = F(z,t) = 1 we can conclude 


c(z, t) < 


cosh7i(l — z) 
cosh(7i) 


(43) 


(44) 


We note that if a(f) is constant in time then 71 defined in (34) is exactly 70 defined in 
(27) and thus F{z.t ) = c(z, 0) proving (35). 

□ 
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With Lemma 5 proven we can find a better lower bound on the pore function t/(z, t ) than 
the one given in Lemma 4. 

Lemma 6 

The pore function rj(z,t) is bounded from below by 

T](z,t) > 7/0-/ ^0{r)F{z,T)dT (45) 

Jo 2 

where F(z,t) is given in (38). 

Proof 

By Lemma 5, c(z, t) < F(z,t). Substituting in (32) we get (45). 


□ 

The above lemmas provide us with an insight into the behavior of the solution. We 
summarize the main points as follows: 

• For any fixed time t, the pore function rj(z,t) is a monotonically increasing function 
(in z), achieving the value 7/(0, t) (defined in (28)) at z = 0. The function r](z,t ) for 
any time is bounded from above by the initial value t/ 0 , and from below by (45). 

• At the time tg the solution becomes discontinuous. In fact integrating (22) between 
z = 0 and z = 1 we get 

-/MO- t)) d - ™J ~\ z =o = J o a 2 (t)rj(^t)c{^t)df (46) 

Since the RHS is positive dc ^ | 2= o must tend to infinity!! 

• The concentration function c(z, t ) is bounded from above by (44) and is monotonically 
decreasing as a function of z for any t. The derivative at z = 0 is becoming more and 
more negative. 

• The function 7 /(z,f) is a decreasing function in time. 

4 Optimization 

After establishing the properties of the solutions, we turn to discuss the concept of a suc- 
cessful process. In CVI, it is generally desirable to produce a solid with a relatively uniform 
porosity distribution. 

The process obviously terminates if 7/(0, t) is small, and we would like rj(z,t) to be uni- 
formly small, i.e. we do not want tj(z, t) to have a large variation in space. Consider therefore 
the average of c 
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Definition: 


m = (47) 

Jo 

By a successful process we mean a process such that e(t) < if for a given target value If 
for some tj < t@. We will show that for each target value If there are functions a(t),0(t) 
such that the process is successful. The question is how to pick a(t),0(t) such that the final 
time tf is minimal. 

We start by observing 

Lemma 7 

Let e(t) defined in (47). Then e(t) is a monotonically decreasing function of t 
Moreover 


de(t) 

dt 


m 

aHt) 




dc(z, t) 

~dz 


z= 0 


(48) 


Proof 

Upon differentiating i(t) we get 


de(t) 

dt 


= 2 / T) (£,t)T)t{£,t) 

J 0 

( 49 ) 

= -M)j 

f r?(£,t)c(£,t)d£ 
0 

(50) 

m 

a 2 {t) 


(51) 

= U= ° 

(52) 

< 0 


(53) 


For the first two steps we used the differential equations (21) and (22), the last step used 
the fact that c — is negative (see Lemma 3 (31)). 


□ 

Lemma 7 implies that the average concentration e is decreasing in time for any parameters 
a(t), d(t). This however does not necessarily yield that the process is successful and that 
there is a time tf <tp in which e < 7/. We need a better estimate of e. This is given in the 
next lemma. 

Lemma 8 

Let e(t) be the concentration average defined in (47). The following estimation holds 
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«(*) < «o - / 
io 


f^\AM7w5^»ta nh a(r) J^^dr 


where r?(0, t) is given in (28). 

Proof 

Let 


d(z,t ) — 


cosh-y^l — z ) 
cos l 172 


g(M) 


Clearly the function d(z,t) satisfies 




Q 2 (t)r?( 0 , i)d(2, i) 


<*(o,t) = i , ^(M) = o ( 5? ) 

We multiply now (22) by d(z,t), (56) by c(z, t), subtract and integrate between 2 = 0 
and 2=1. Taking the boundary conditions for c(z,t), d(z,t) into consideration we get: 


-/OKo,0)^U + /MM))^U. = 


a 2 (t) J o - 77(0 ,t))c(t,t)d{t,t)d£ 

+ [ urn, t) - mo, d) ;) m . | t} . 


The right hand side of (58) is positive leading to 


dc(z,t) t ^ dd(z,t) l 
dz iz=0 - d~z l,=0 


dz ~ dz |2 ~ u 

= -7 2 tanfi(7 2 ) 

We turn now to (48) and integrate (in respect to the time t) to get 




dc{z,T ) 


Substituting (59) into (61) yields (54). 
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□ 

We will show now that for any target value l s there are functions a{t ) and /3(t) such that 
the process is successful. In fact we will show that a successful process can be obtained by 
choosing time independent a(t) = a and 0(t) = 0. In this case (28) reads 

v(0,t) = Vo-^0t (62) 

and the terminal time is given by 

t 0 = 

In the following Lemma we show that we 
successful process. 

Lemma 9 

Let rj(z, t) be the solution of (21), (22) with fixed a and (3. Let l be defined in (47). Then 

= 0 (64) 


2t?o 

0 


(63) 


can always choose a small enough to have a 


Proof 

We start by using (54) for fixed a and 0 and the expression (62) for rj( 0, t) to get 


-, 2 Vo, ^ a fn , tanh(g{T,a )) 

-( T ) < <0 -J 0 M0,r) dr 


Where 


g(t,a) = a (t), 

We introduce now a change of variables 


[ 7(0,*) 

^ /(»?(0, <)) 


£ = 


0T 

2% 


To get 


1^) < to - 2vl [\l - 0 

p Jo 


tanh(G(^,a)) 

G(S,a) 


dfi 


(65) 


( 66 ) 


(67) 


We use now the Lebesgue theorem and the fact that (?(£, a) tends to zero as a tend to 
zero to conclude that 


lim Q -> 0 e 


^) < eo-2e 0 [\l-S)dt = 0 
p Jo 


( 68 ) 


This proves the Lemma. 
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□ 


The parameters a and (3 are not independent. In fact (3 is proportional to a - 2 see (16). 
Thus if we choose these parameters as constants as in Lemma 9. the time needed for the 
completion of a successful process is very long. We would like to minimize that time. We 
therefore formulate the following minimization problem. 

Problem: 

Let r](z,t),c(z,t) be the solution of (21-25) for some a(t), 3{t). Let the average concen- 
tration be defined in (47). Given a target concentration value If we define the final time 
tf(a,0) as the first time such that 

*(*/) < e 7 

find min Q(t) , 5 ( t) t f (a,{3). 


We did not solve this optimization problem, however Lemma 9 provides us with some 
insight to the minimization process. In fact it is easier to treat the upper bound on the 
average concentration given in (54). 

In the following we will show an example in which some savings in the time needed to 
reach the target value e(tf) is realized by choosing a and 0 that are time dependent. 

Consider the special case 

3{t) = 2 a 2 {t) 

and 

/(*?) = *?' 3 

This expression for f(rj) corresponds to the case were N k « 1 and 9 = 77 -1 . Physically, 
the small value of N k means that Knudsen diffusion contributions are not significant. As the 
porosity is filled by the matrix phase this assumption will invariably break down, however, 
with relatively large pressure and/or pore sizes, N k will only become significant at the end 
of the process. In this case, the analysis presented below (i.e., where N k was ignored) will 
be reasonably accurate. 

In this case we may rewrite the estimate (54) for the average concentration in the form 


where 


m 


< 




_ r B{TH0 . T) ^±a dT 

Jo 72 


72 (<) = 


* 7 ( 0 , t) 


(69) 
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and the pore function 77 ( 0 , t) is given in Lemma 2. From this Lemma we can deduce that 

8 = -2^(0, t) (70) 

* 

We will look now for a design function a(t ) such that 72 is independent of time. 

This leads to an explicit expression for a(t), in fact solving 

« 2 (* ) 02 

vM 2 


and taking into account that 

a 2 (t) 



we get 

77 ( 0 , t) 

Vo 

1 + r) 0 R 2 t 

(71) 


a{t) 

\R\vo 

1 + rjoR 2 t 

(72) 


We can now substitute the result in (69) to get 


< Vo ~ 


tank 72 
72 


\M 

(1 + T] 0 R 2 t) 2 



(73) 


The parameter R can be chosen to minimize the time needed for e(t) to be less than the 
target value. This approach leads to some savings, as described below. 


5 Numerical Results 

We have solved numerically the following set of equations 

(rfc z ) + a 2 r]c = 0 

Vt = -8c 


subject to the initial condition 
and boundary conditions 


77 ( 2 , 0 ) = 1 


c( 0 , t) = 1 c 2 (l, t) = 0 


(74) 

(75) 

(76) 

(77) 


To discretize (74) we used the pseudospectral Chebyshev method. In this method we choose 
the grid points to be 


Xj — 


1 + cos(?4) 

7— 1 0 <J < ft 


(78) 


14 



N is the total number of grid points. 

The spectral differentiation matrix takes the value of a given function at the gridpoints 
Xj and yields the values of the derivative of the interpolation polynomial at these points. 
The points Xj are the nodes of the Gauss Lobatto Chebyshev quaderature formula. The 
matrix can be written explicitly: 


Djk = 

Djj = 
Aw = 


1 c, (-!)>+* 

2 c k sin sin ^(-j+fc) 


-1—3 

2 sin^j) 

—Dnn = 


2A r2 +l 

6 


3*0, N, 


(79) 


We apply the matrix Djk twice, once for the vector c taking into account the boundary 
condition c(0, t) — 1 and then to rfc z now taking into account that c z (l,t) = 0. This yields 
a linear system for the values of c(xj, t). 

In the next stage we update 77 by the standard fourth order Runge-Kutta scheme. 


• We first ran the scheme with constant a ( (3 = a 2 ) with the requirements that the 
process stops when e is less than 5% of the original value, and also that e(l,t) will 
be less than 8%. The largest a (and therefore the shortest time to reach the end of 
the process) that satisfied the conditions was a = .235. The total time needed for the 
process to settle was 

_ 2.458 

tfl ~ W) 


• We then used formula (72) to choose a(t) where c*(0) = .3. Again / 3(t ) = a 2 (t). The 
process terminated successful}' at time 


tf2 — .963t/i 


• A better strategy was to use (72) until a became smaller than .25 and then to raise it 
to .27 and reapply (72). Here too /3(t) = o 2 (t). The process terminated here at 

t /3 = .935t/i 

Whereas the savings do not seem extremely significant, they indicate that one can get better 
results by varying the design parameters a and (3 with time. 

• Another possibility is to change the relationship between a and (3 in time. This is 
motivated by Lemma 9. It is clear from this Lemma that for the process to converge 
we need a small. It is also clear that the time for the process to settle is proportional 
to (/3) _1 , thus one wants to minimize a to get a successful process and to maximize !3 
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to reduce time. We ran our problem with a(0) = .4 and /?(0) = a 2 (0). We then varied 
a in time according to (72) while holding /? fixed. We got convergence at 


t f 4 = .397*,!!!! 


(the final ct was .2). This result indicates that great savings can be obtained by choosing 
appropriate values for the design functions a(t) and j3(t). 


6 Conclusions 

In summary, the mathematical results presented here are potentially important in several 
ways. First they provide detailed information on the solutions in space and in time, as 
well as upper and lower bounds. The formulation and general analysis provide a basis for 
understanding the infiltration process in terms of only two parameters: a and /?. The 
treatment of the optimization problem makes it possible to assess optimal (i.e. minimal) 
infiltration times concisely. The observation that the minimum time can be dramatically 
decreased by varying a with time is particularly important because of the long times that 
are typically associated with CVI. Varying a while fi is constant requires varying one or 
more of the relevant processing conditions (e.g., T, Cr, or P) in a controlled way. A detailed 
analysis of these possibilities will be presented elsewhere ([12]). 
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